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Abstract. We present the theoretical background and detailed equations for the syn- 
chrotron emission of a shock wave propagating in a relativistic jet. We then show how the 
evolution of an outburst in this shock-in-jet scenario can be analytically described and pa- 
rameterized to be fitted to multi-frequency lightcurves of galactic and extragalactic sources. 
This is done here for the first time with a completely physical description of the jet and the 
shocked gas, while previous studies used a more phenomenological approach based on the 
observed properties of the outbursts. Another interesting addition to previous work is the 
introduction of a low-energy cut of the electron energy distribution that allows for much 
more diverse synchrotron spectral shapes. To demonstrate and illustrate the new method- 
ology, we present results of infrared-to-radio lightcurve fitting of a succession of outbursts 
observed in 1994 in the microquasar Cyg X-3. We find that the diversity of outbursts in 
shape, amplitude, frequency range and timescale can be fairly described by varying only 
the strength of the shock and its build-up distance from the apex of the jet. A rapid build- 
up results in high-frequency outbursts evolving on short timescales, while slowly evolving, 
low-frequency outbursts form and evolve further out in the jet. We conclude by outlining 
future developments, in particular the inclusion of the associated synchrotron self-Compton 
emission at X-rays and gamma-rays. 

Key words, radiation mechanisms: non-thermal - galaxies: active - galaxies: jets - radio 
continuum: galaxies - radio continuum: stars 



1. Introduction 



25 years ago, iMarscher & Geart dl985l) in- 
troduced a shock-in-jet model to explain a 
giant flar e in the quasar 3C 273 observed 
in 1983 dRobson et al.l Il983l) . The observa- 
tions showed a rise of the flux density of 
the synchrotron self-absorption turnover to- 
gether with a slight move of the turnover 
from high- to low-frequencies. This behavior 
wa s not compat i ble wi th the plasmon model 
of Ivan der Laa ij| {[966) describing the adia- 
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batic expansion of a spheroidal blob of plasma. 
It was, however, compatible with a shock 
model taking into account Compton cooling 
of the el ectrons in the initial s tage of its 
evolution (IMarscher & Gear| [l985). Until now, 
this model could not be disproved and was 
found to give a good description of the ob- 
servations of both galactic and extragalactic 
sourc e s of r elativistic iet s (3C 273:lTurler et al J 
I1999L l2000h (3 C 279: ILindfors et al.l l2006h 
(GRS 1915+105:lTurier et al.l2004 (Cyg X-3 : 
ILindfors et alJl2007t iMiller- Jones et al.ll2009b . 

Based on the experience gained on the 
modeling of the flaring behavior of these 
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sources, we describe here a fully physical ap- 
proach that we intent to use in the future. It 
has the advantage to fit directly the physi- 
cal jet properties (electron energy distribution, 
magnetic field, etc.) rather than the observ- 
ables (frequency and flux density at the spec- 
tral turnover, etc.). This change will also al- 
low us to soon incorporate the associated syn- 
chrotron self-Compton component in a self- 
consistent way. We present the main equa- 
tions of the shock model starting from stan- 
dard synchrotron theory in Sect. [2] and the 
method of shock-in-jet modeling in Sect. [3] We 
then present results obtained for Cyg X-3 in 
Sect. [4] and conclude with future perspectives 
in Sect. |5] 

2. Theory 

Synch rotron theory as derived by Pach olczvkl 
(1970), for instance, is relatively com- 
plex. Let's assume a homogeneous syn- 
chrotron source with a powerlaw electron 
energy/impulsion distribution given by 
n(y) = Ky~ p for y min < y < y max , where 
n is the electron number density in cirT 3 
and y is the Lorentz factor of the relativistic 
electrons, or more precisely the product /3y 
for mildly relativistic electrons. This results in 
equations for the emission s v and absorption 
k v coefficients that have their simplest form 
when expressed as a function of the cyclotron 
frequency v# = eB/(2nmc) where e and m are 
the charge and mass of the electron and B is 
the component orthogonal to the line of sight 
of an uniform magnetic field in the source. 
The equations are: 
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where g s (p) and g K {p) are a product of Euler 
gamma functions, T, given by: 
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The latter functions are slowing varying in the 
range of interest and have both a minimum at 
g s (p* 2.6) ~ 8.2 and at g K {p* 0.9) « 26.4. 

Solving the differential equation of the ra- 
diative transfer then gives the specific intensity 
I Y in erg s _1 cirT 2 Hz -1 sr~' as: 

/v^d-e-^^^a-e-). (5) 
Kv 8k(P) v^ 2 

One recognizes the typical synchrotron spec- 
trum with an optically thick part (t v » 1) go- 
ing as v 5//2 and the optically thin part (r v <K 1) 
following a v^ p ^ 1 ^ 2 dependence since /' hln = 
xs v , which is obtained by replacing (1 - e~ Tr ) 
by t v — x k v in Eq. (0, where x is the thickness 
of the source along the line of sight. We note 
that this standard synchrotron spectrum strictly 
holds only if there is no significant radiative 
cooling of the electrons and if v m [ n tends to 
and v max tends to infinity, where v min = v B y 2 min 
and v max = VBT^ax- A synchrotron spectrum can 
have a much different shape if radiative cooling 
is important or if v m ; n happens to be higher than 
the frequency of synchr otron self-absorption , 
v a bs, as illustrated by e.g. lGranot & Saril(l2002i 
Fig. 1). In particular, in the slow cooling case, 
the optically thin spectral slope will break to 
a v 1 / 3 dependence at frequencies between v a b s 
and Vmin- 

For the emission of a shock wave in a 
quasar, we consider a source located at a lu- 
minosity distance Di - corresponding to a red- 
shift of z - with a slab geometry of thickness x 
and radius R, corresponding to the half-width 
of the jet (see Fig. [TJ. The source is mov- 
ing with a speed Be corresponding to a bulk 
Lorentz factor of F = (1 -B 2 y 1 ^ 2 in a direction 
making a small angle 9 with the line-of-sight, 
which results in a Doppler boosting factor 6 — 
r _1 (l - B cos 9y l . The optically thin flux den- 
sity F v hm = O /' hln = Q. x s v of a uniform source 
sustaining a solid angle Q. = jtR 2 (1 +z) a ID 2 l is 
then given by: 



( j(p+5)/2 ^2 
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where x, R, K and B are expressed in the frame 
of the moving source, whereas v and F v are the 
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adiabatic synchrotron Compton 
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Fig. 1. Schematic representation of a propagating shock wave in a r elativistic jet and the three -stage evolu- 
tion of the associated synchrotron outburst acco rding to the mode l of Marscher & Gear (1985) and with the 
modification of the Compton stage proposed by Biorn sson & Asla ksen (2000) (dashed line). Here s = p. 



observed quantities. We can also derive the fre- 
quency of synchrotron self-absorption by set- 
ting r v (v abs ) = x abs Ky(v a bs) = 1, which gives: 



6 / e g K (p) r (p+Xi/z\ 



v 2/0+4) 



1 +z \ 16rac 



(7) 



The two last equations expressed as pro- 
portionaliti es are the starting poin t of the shock 
model of lMarscher & GearlTl985h . The critical 
parameter that defines the well-known three- 
stage evolution of the outburst (see Fig. []]) is 
the thickness x of the emission region along 
the line of sight. The evolution of a spher- 
ical plasmon expanding adiabatically in 3- 
dimensions is obtained by replacing x by R 
and by let ting K evolve with R'^ 2 ^ and B 
with R~ 2 ( Ivan der Laanl Il966h . The result is 
a decrease of flux density, while the spectral 
turnover is moving towa rds lower frequencies. 
Marsc her & Gear! (1 19851) assume the source to 
be a cylinder of thickness xr - faR - a con- 
stant fraction of R - in a conical jet ex- 
panding in the two directions perpendicular to 
the jet axis such that K oc /r 2( ' J+2)/3 and with 
the magnetic field to be preferentially orthogo- 
nal to the jet axis (B oc R l ) which leads to a 
slightly shallower decrease of the flux density, 
but the trend is the same. To explain the rise in 
flux density observed in the early phases of the 
the 19 83 outburst in 3C 273. iMarscher & Gearl 
d 19851) had the idea to strongly limit the source 



volume at the onset of the outburst. This is 
achieved by considering that x <K xr at the 
beginning of the shock evolution because the 
thickness is limited by the distance electrons 
can travel from the shock front before they cool 
down throu gh radiative losses, such that x (v) = 
ffmicfrnniM. lBiornsson & Aslaksenl d2000l) pro- 
po se a modification of the e xpression of x used 
bv lMarscher & Gearl d 19851) as: 



x(v) 



rel 



lcr T U, 



1/2 

B + U S \l+Z> 



1/2 



,,-1/2 



(8) 



where <r T is the Thomson cross-section and 
B K \ is the average speed - in units of c - of 
the electrons relative to the shock front, which 
takes a value o f 0.3 for a relativistic shock 
dBiornssonll2010lPl The dominant cooling can 
either be inverse-Compton on the virtual pho- 
tons of the magnetic field with an energy den- 
sity Ub = B 2 l(%n) (synchrotron cooling) or 
on the actual synchrotron photon with an en- 
ergy density of Us — (4-n/c) J y I v dv (first-order 
Compton cooling). In any case, this radiative 
cooling results in a steeper optically thin slope 



1 The factor 3/8 in Eq. l[8j corresponds roughly to 
the cooling time needed for the final (cooled) emis- 
sion frequency to be at half the initial frequency, 
Vf/Vi = (7f/7i) 2 = 1/2, rather than yt = yJ2, which 
corresponds to a factor 3/4 (Bjornsson C.-L, private 
comm.). 
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following a v ~ p ' 2 dependence compared to the 
adiabatic cooling, as can be seen by inserting 
Eq. ([8]) into Eq. ©. But even in the adiabatic 
stage, a cooling break remains at v coo i - defined 
by x(v CO oi) = Xr - resulting in a spectral index 
of -p/2 at higher frequencies. 

The next step is to insert Eq. (O for jt a b s = 
xivabs) into Eq. ^} to get: 



/ (p+5)/2 _ 3e 2 cB xA g K (p) 
Vabs ~ 128cr T (t/B + f/ S ) 



Kv 1 



,(p+3)/2 



(9) 



where we introduce V — v ( 1 +z) 15 to denote the 
frequency measured in the frame of the moving 
source. 

The main difficulty comes now with the 
evaluation of Us- The main contribution to Us 
comes from the optically thin part of the emis- 
sion I v hm = x(v) s v , so that we can calculate Us 
by integrating the optically thin spectrum ex- 
tending between vi and V2 by using Eqs. ([T} 
and ©, as: 

uI + UbUs = ^Hl^-Kyf f W(10) 
For p + 2, the integral in Eq. (ITOb is equal to 



2(v: 



/(2-p)/2 ,/(2-p)/2 



-V 



) / (2-p) and then the second 



order equation can be solved for Us- While V2 
is always v max , Vi is the highest of the two fre- 
quencies v a bs or v m i n . Whether this is the case 
is however not known a priori since v a bs de- 
pends on Us in Eq. (|9j. The way to proceed 
is to calculate Us with v\ set to v ra j n to get a 
first guess of v a b s , which we call v a b s ,i. Even 
if it turns out that v a b Sj i > Vmin, we can still 
stay with this result for Us in most cases. It is 
only if in addition p > 2 and Us is compara- 
ble or greater than Ub that taking v a b s , instead 
of Vmin, for v\ would have a significant effect. 
In the latter case, the following expression de- 
rived from Eqs. (|9]l and (TT~Qb was found to be a 
good approximation for the right-hand side of 
Eq. dlOll when p > 2: 



3mne 2 p K i K v 2 B 
8(p-2)cr T g-Hp) 
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Once Us is derived, it is simple to calcu- 
late ;t a bs from a relation obtained by inserting 
Eq. (0 into Eq. ([8]). We then limit x a b s to its 



maximum value of xr applying to the final de- 
cay stage and calculate v a b s through Eq. and 
^ abs = n hin (v abs ) through Eq. ©. 

3. Method 

Applying the theory outlined in Sect. [2] to 
multi-wavelength observations of quasars and 
microquasars is the next challenge. For this, 
we need to define the temporal evolution of the 
physical parameters K and B; of the spectral 
turnover defined by v a b s and F a b s ; and of the 
frequency of additional spectral breaks at v m i n , 
and v coo i- 

We consider a jet that is not necessarily 
conical, but with an opening radius R increas- 
ing as a powerlaw of the distance along the jet 
X as R = Ro(X/X(,y, where Xq is an arbitrary 
normalization length. If r < 1 we get a colli- 
mating jet, whereas if r > 1 we get a trumpet- 
shaped jet with decreasing collimation. 

We can then assume powerlaw depen- 
dences of the electron number density K - 
Ko(R/Ro)~ k and of the magnetic field B = 
Bq(R/Ro) along the undisturbed jet. For an 
adiabatic jet flow expanding in two dimensions 
perpendicular to the jet axis, we have A: at j = 
2(p + 2)/3, which can be seen as a lower-limit 
if there are also significant radiative losses. The 
value of b shall be between 1 and 2 correspond- 
ing, respectively, to the perpendicular, B ± , and 
parallel, By, components of the field. There are 
however two arguments to prefer a value of 
b = 4/3. First by assuming a turbulent mag- 
netic field one gets B ± = 2B\\ - rather than 
B± = B\\ - that corresponds to b — 4/3, sim- 
ply because there are two dimensions of space 
perpendicular to the jet axis and only one par- 
allel to it. Secondly, by assuming equipartition 
between the energy densities of the electrons, 
U e - mc 2 J n(y)ydy, and the magnetic field, 
U B = B 2 /(8tt), one gets B 2 /K oc rUp-VI* and 
thus b ecj = (k/2)-(p-2)/3, which is equal to 4/3 
for k = k n< i- 

We note that in t he sh ock-in- jet scenario 
of IMarscher & Gearl d 19851) . the electrons are 
already accelerated in the undisturbed jet and 
the shock wave mainly compresses the flow 
locally without strongly accelerating particles. 
The strength of the compression factor, r\ - 
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defined as the ratio of electron number den- 
sities in front of and behind the shock front 

- is a good candidate to account for differ- 
ences from one outburst to the next. We there- 
fore leave r\\ as a free parameter for each out- 
burst i. The evolution of B and K for out- 
burst i is thus given by B, = rjiBo(R/Ro) and 
Kj = T]jKo(R/Ro)~ k . The energy gain of each 
electron as it crosses the shock front, has a 
minimal value of £,d = rf^ in the case of an adi- 
abatic acceleration process ( Mars cher & Gear! 
1 1985b . By taking this as the baseline, we get 
y m = 7 1 1 i /3 j mfi (R/R )- (k - 2mp - 1 \ where the sub- 
script 'm' stands for 'min' or 'max' and where 
the exponent of R - equal to -2/3 for k = £ a d 

- is deduced for a two-dimensional expansion, 
i.e. with the electron number density decreas- 
ing as n — J riydy oc R 2 . This leads to v m = 
(6/(l+z)) VBylocrf^R-b-W-W-P-^, where we 
assumed a constant Doppler factor 6. The evo- 
lution with R of the two final spectral breaks, 
v C ooi and v a bs, is more complex as they depend 
on Us, which can be obtained from Eq. ( TTOb . 

As we are interested in the evolution of an 
outburst with time rather than with R or X we 
define an origin of time for the jet flow at the 
apex of the jet, tj et (X — 0) = 0. In the frame 
of the black hole, we have X - Be f™, which 
gives X = BcTd fj e t/(l +z), by expressing the 
time in the observer's frame. To account for a 
possible acceleration or deceleration of the jet 
flow, the simplest is to express it as a power- 
law dependence on X, as BY — BoTo(X/Xo) g 
and 5 = 5o(X/Xo) d , where we linked B and F 
to account for both mildly and strongly rela- 
tivistic jets and we separated the Doppler factor 
to further account for the emission of a curved 
jet. For a straight jet, we have qualitatively that 
d x for B F < 1 , whereas d depends on g for 
BT>\ as d ~ g for 6 < l/F and as d » —g for 
> l/F. We can then integrate the expression 
of dtjet = (1 +z)/(cBF6)dX from zero to fj et to 
get: 



R=R o[§- \= R " 



'jet 
fjet.O 



r/(l- g -d) 



(12) 



where f jet;0 = (1 +z)X /(cB Fo6o(l-g-d))- This 
relation links R to the observable fj et = t - 
tx=o, where tx=o can be fitted for each outburst 



present in the observed lightcurves. It is how- 
ever more natural to fit the time of the peak 
of the outburst, t p j, and the corresponding dis- 
tance along the jet, X p j, for each flare i and to 
use this to calculate tx=o,i- The choice of X p as 
a free parameter for any outburst in addition 
to rji is motivated by the finding that the typ- 
ical distance where the shock evolves seems 
to be the main driver of the different proper- 
ties of individual outbur sts found in Cyg X-3 
( Mille r- Jones et a 1. 2009). As it is not physical 
to set the compression rj to zero until reach- 
ing X p j and then abruptly let it jump to its final 
value T]i, we assume a linear increase of rj with 
X from the apex of the jet to X p j. 

To summarize, the parameters we fit to a 
set of multi-frequency lightcurves are some of 
the indices p, r, k, b, g, and d depending on the 
assumptions, as well as some of the normaliza- 
tions K , B ,B Fo, Tmin.o, Tmax.o, and the factor 
/«. With this and the fixed value of Xq, we can 
get /?o -X{) tan <po based on the jet o pening half- 
angle 0o * 10° /To deduced by iJorstad et al.l 
(2005). 5a also derives from To by using the jet 
angle to the line-of-sight, 9, taken from stud- 
ies of superluminal motion in the considered 
sources. We relate the observed time in the 
lightcurve t to the radius of the jet /? through 
Eq. (fT2l and we can thus calculate the evolu- 
tion with observed time of K, B and x a b s and 
thus of the synchrotron spectrum completely 
defined by v abs , F abs , and Vcooi- 

Except for the high-energy end of the spec- 
trum - which we impose to cut-off exponen- 
tially at v max - the modeling o f the shape of 
the spe ctral breaks is taken from lGranot & SarH 
(2002, Eq. (4)). We made, however, the simpli- 
fication of setting the sharpness of the break 
to s = 3/(B\ -Bi), where the factor 3 has been 
chosen to reproduce well the synchrotron self- 
absorption turnover, which is the most impor- 
tant break. Although this is likely less accu- 
rate physically, taking an inverse dependence 
of s on (B1-B2) has several advantages. Firstly, 
it has additive proprieties such that twice the 
same break at the same frequency is equivalent 
to a twice sharper break of 2(J3\ —fh), thus en- 
suring smooth transitions when a break crosses 
the frequency of another break. Secondly, it 
also works well for negative (concave) breaks 
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and a break of ySi —fh = is equivalent to no 
break, whereas otherwise it changes the nor- 
malization of the powerlaw. 

The model lightcurve of an outburst is con- 
structed from its evolving spectrum by extract- 
ing values at different times, but at a given fre- 
quency. The fitting of a dataset is done simul- 
taneously for all available lightcurves, but only 
on a subset of all free model parameters at a 
time. A series of several fits covering all the 
parameter set, is repeated iteratively until con- 
vergence. The number of outbursts and their 
approximate onset time and amplitude is often 
defined manually beforehand to guide the itera- 
tive fitting process. We implemented the possi- 
bility to fit the data with additional constraints 
to favor solutions with parameter values closer 
to expectations, e.g. for a simple conical (r = 1) 
and adiabatic (fe a( j) jet flow with equipartition 
(b eq ) and/or with a minimal dispersion among 
the specificities X p j and 77, of the different out- 
bursts i. Moderate random variations of the un- 
certainties associated to each data point in a 
lightcurve before consecutive fit iterations was 
found to be a good trick to ease convergence in 
blocked situations, i.e. when the fit is stacked 
in a local minimum in x 2 , but far from the ab- 
solute minimum. 

Finally, we note that the contribution from 
the underlying, undisturbed jet flow can be de- 
rived from the physical parameters of the jet 
by setting X p = and r\ = 1, while the cumu- 
lative contribution of decaying outbursts peak- 
ing before the start of the considered dataset 
can be obtained by considering a series of out- 
bursts with average values of X p and 77, spaced 
by the typical time interval between consecu- 
tive events. 

4. Results 

As an illustration of the physical modeling de- 
scribed above, we present here a fit to the 
dataset of C yg X-3 observed in February- 
March 1994 ((Fender et all 119971) . We chose 
this dataset because it includes a wide va- 
riety of outbursts differing in amplitude, 
time- and frequency-scale, as derived by the 
more phenomenolo gical approach used by 
iLindfors et all (120071) . The resulting fit to the 
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Fig. 3. Evolution of the average model outburst, a) 
Evolution with time of the peak (red line with ar- 
rows) of the synchrotron spectum (thin gray lines) 
resulting in a peak flux reached at different frequen- 
cies (green line). The peak of the overall outburst's 
evolution is shown by a star symbol, while circles 
show the position of this peak for each of the 14 in- 
dividual outbursts, with same colors as in Fig.[2] The 
four vertical lines at the top indicate the frequencies 
of the lightcurves in the dataset. b) An areal view 
on the three-dimensional model outburst in the flux 
versus time and frequency space. Thin gray lines 
are contours of equal flux density, whereas the other 
lines and points are as in the upper panel. 



three radio-band lightcurves and the few data 
points in the infrared K-band is shown in 
Fig. [2] We obtain a fair description of the 
dataset with a set of 14 self-similar outbursts, 
and thus demonstrate that our assumptions of 
varying only the strength of the shock and its 
build-up distance is enough to describe the ob- 
served differences among the outbursts. The 
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Fig. 2. Fit of a series of 14 model outbursts - shown by different colors - to the flaring episode observed in 
Cyg X-3 during February-March 1994. 



colored points in the lower panel of Fig. [3] 
show that the trend from short-lived, high- 
frequency peaking outbursts to long-lasting, 
low-frequency peaking outbursts is well repro- 
duced by varying the distance, X p , along the jet 
where the shock compression reaches its max- 
imum. We find an average distance among the 
outbursts of Xp = 1.8 X 10 15 cm with typical 
fluctuations by a factor of 2.5 from one out- 
burst to the next. 

We find here that the assumptions of a con- 
stant jet flow (g — d — 0) expanding adiabati- 
cally (£ ac j) in a conical (r = 1) jet, and with the 
preferred index b eq = 4/3 for the decrease of 
the magnetic field was fine for the considered 
dataset. Because of degeneracy among some of 
the parameters, we further fixed the value of 
to 10" 2 ,/?oro to 1 typical for a moderately rela- 



tivistic jet with an ass umed jet angle to the line - 
of-sight of 6= 10.5° dMiller- Jones et al.ll2004l) . 
and impose an average comp ression factor of 
p = 5 , according to r\ £ 6 (Mars cher & Gearl 
1985). Therefore, the only jet properties that 
are free to vary are the index, p, of the elec- 
tron energy distribution and the normalizations 
at Xq = 10 16 cm of its low-energy bound, ymin,o> 
as well as the corresponding normalizations of 
the magnetic field Bo and of the electron den- 
sity K(). The exact value of y maXi o was found to 
be irrelevant and we thus fixed it to 10 5 . 

The best fit parameters are an electron in- 
dex of p = 2.01 and a low-energy limit of 
Tmin.o =4.0 for the electron energy distribution, 
with a number density of Kq = 3.4 cm -3 , and 
a magnetic field of Bq - 54 mG. These quan- 
tities apply to the underlying jet flow at a dis- 
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tance of Xq = 10 cm from the apex of the jet. 
For an average outburst peaking at a derived 
distance of X p = 1.8 x 10 15 cm from the apex 
of the jet, we get peak values of y m in, P =21, 
K p = 1.7 X 10 3 cm 3 and B p = 2.5 G. The cor- 
responding energy densities of electrons and 
magnetic field are I/ e ,p = 1,2 X 10~ 2 ergcirT 3 
and £/B,p = 0.25ergcirT 3 , respectively. 

The upper panel of Fig. [3] shows the evo- 
lution with time of the spectrum of the aver- 
age outburst. Although this looks very simi- 
lar to the typical evol u tion o f the shock model 
of iMarscher & Gear] (1 19851) without a syn- 
chrotron plateau stage (see Fig. Q]), the initial 
rise of the spectrum turnover is of a completely 
different nature. The parameters we derive do 
not lead to significant radiative losses (we al- 
ways have v coo i » v a b s ), so that we are in the 
adiabatic cooling stage throughout the outburst 
evolution. It is mainly the increasing compres- 
sion j] and secondly the presence of the low- 
energy break at v m i„ that define, actually, the 
initial rise of the spectrum. 

5. Conclusion 

We presented a new, physical modeling and 
parametrization of outbursts in relativistic jets. 
We applied the model to a rich dataset of Cyg 
X-3 and derived the physical conditions in the 
jet. We found that standard assumptions of a 
conical, adiabatic jet flow seem adequate and 
that conditions are such that radiative cooling 
is negligible. The next step is to fit the model 
parameters to other datasets as it has been done 
in the past with a more phenomenological ap- 
proach. We will then be able to compare phys- 
ical conditions in the jets of different sources 
from microquasars to blazars. 

An interesting addition to the modeling 
would be the inclusion of the associated syn- 
chrotron self-Compton emission in the X- and 
gamma-ray spectral domain, where Swift/BAfT 
and Fermi are currently providing an almost 
continuous monitoring of bright blazars. This 
has already been done in a simplified way 
by iLindfors et alJ (|2005) and would be par- 
ticularly interesting for sources where multi- 
ple inverse-Compton scattering might play an 
important role. We would then build on the 



theoretical study of Bjornsson] d2010h to con- 
struct a self-consistent description of the multi- 
frequency emission of shock waves in relativis- 
tic jets from the radio to the gamma-rays. 
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